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Abstract 

We investigate the properties of the hadron-quark mixed phase in compact stars using a Brueckner- 
Hartree-Fock framework for hadronic matter and the MIT bag model for quark matter. We find 
that the equation of state of the mixed phase is similar to that given by the Maxwell construc- 
tion. The composition of the mixed phase, however, is very different from that of the Maxwell 
construction; in particular, hyperons are completely suppressed. 
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1. Introduction 

It is well known that hyperons appear at several times normal nuclear density and lead 
to a strong softening of the equation of state (EOS) with a consequent substantial reduc- 
tion of the maximum neutron star mass. Actually the microscopic Brueckner-Hartree- 
Fock approach gives much lower masses than current observation values of ^ 1.5Mq. 

On the other hand, the hadron-quark deconfinement transition is believed to occur in 
hot and/or high-density matter. Taking EOS of quark matter within the MIT bag model, 
the maximum mass can increase to the Chandrasekhar limit once the deconfinement 
transition occurs in hyperon matter P][2] . Since the deconfinement transition from hadron 
to quark phase may occur as a first-order phase transition, the hadron-quark mixed phase 
should appear, where charge density as well as baryon number density is no more uniform. 
Owing to the interplay of the Coulomb interaction and the surface tension, the mixed 
phase can have exotic shapes called pasta structures . 

The bulk Gibbs calculation of the mixed phase, without the effects of the Coulomb 
interaction and surface tension, leads to a broad region of the mixed phase (MP) [4]. 
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However, if one takes into account the geometrical structures in the mixed phase and 
apphes the Gibbs conditions, one may find that MP is considerably limited and thereby 
EOS approaches to the one given by the Maxwell construction (MC) |3j. 

In this report we explore the EOS and the structure of the mixed phase during the 
hyperon-quark transition, properly taking account of the Gibbs conditions. 



2. Numerical Calculation 



The numerical procedure to determine the EOS and the geometrical structure of the 
MP is similar to that explained in detail in Ref. [3]. We employ a Wigner-Seitz ap- 
proximation in which the whole space is divided into equivalent Wigner-Seitz cells with 
a given geometrical symmetry, sphere for three dimension (3D), cylinder for 2D, and 
slab for ID. A lump portion made of one phase is embedded in the other phase and 
thus the quark and hadron phases are separated in each cell. A sharp boundary is as- 
sumed between the two phases and the surface energy is taken into account in terms of 
a surface-tension parameter a. The energy density of the mixed phase is thus written as 



1 



aS 



where the 



volume of the Wigner-Seitz cell Vw is the sum of those of hadron and quark phases Vh 
and Vq, S the quark-hadron interface area, en, eg and eg are energy densities of hadrons, 
quarks and electrons, which are r-dependent since they are functions of local densities 
/9a(r) (a — n,p, A,'E~ ,u,d, s, e). The Coulomb potential Vc is obtained by solving the 
Poisson equation. For a given density ps, the optimum dimensionality of the cell, the 
cell size Rwj the lump size R, and the density profile of each component are searched for 
to give the minimum energy density. We employ a — 40 MeV/fm^ in the present study. 

To calculate ch in the hadron phase, we use the Thomas-Fermi approximation for the 
kinetic energy density. The potential-energy density is calculated by the nonrelativistic 
BHF approach [T] based on microscopic NN and NY potentials that are fitted to scattering 
phase shifts. Nucleonic three-body forces are included in order to (slightly) shift the 
saturation point of purely nucleonic matter to the empirical value. 

For the quark phase, we use the MIT bag model with massless u and d quarks and 
massive s quark with = 150 MeV. The energy density eg consists of the kinetic 
term by the Thomas Fermi approximation, the leading-order one-gluon-exchange term 
[5] proportional to the QCD fine structure constant as, and the bag constant B. We here 
use B = 100 MeV/fm'^ and as = to get the quark EOS which crosses the hadronic one 
at an appropriate baryon density. 



3. Hadron-Quark Mixed Phase 



Figure [T] illustrates an example of the density profile in a 3D cell. One can see the non- 
uniform density distribution of each particle species together with the finite Coulomb 
potential; charged particle distributions are rearranged to screen the Coulomb potential. 
The quark phase is negatively charged, so that d and s quarks are repelled to the phase 
boundary, while u quarks gather at the center. The protons in the hadron phase are 
attracted by the negatively charged quark phase, while the electrons are repelled. 

Figure [2] compares the resulting EOS with that of the pure hadron and quark phases. 
The thick black curve indicates the case of the MC, while the colored line indicates the 
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Fig. 1. Density profiles and Coulomb potential Vc 
within a 3D (quark droplet) Wigner-Seitz cell of 
the MP at pB = 0.4 fm~^. The cell radius and the 
droplet radius are R\y = 26.7 fm and R = 17.3 
fm, respectively. 



Fig. 2. EOS of the MP (thick curves) in com- 
parison with pure hadron and quark phases (thin 
curves), hadron (ps < 0.44 fm~^) or quark 
{pB > 0.44 fm~^) phases. Each segment of the 
MP is chosen by minimizing the energy. 
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Fig. 3. Particle fractions in the MP by the full calculation (left panel) and the MC (right panel). 

MP in its various geometric realizations starting with a quark droplet structure and 
ending with a bubble structure. Note that the charge screening effect, combined with the 
surface tension, makes the non-uniform structures mechanically less stable and limits the 
density region of the MP JJJ . Consequently the energy of the MP is only slightly lower 
than that of the MC. However, the structure and the composition of the MP are very 
different from those of the MC, which is demonstrated in Fig. [21 where we compare the 
particle fractions as a function of baryon density in the full calculation (left panel) and 
the MC (right panel). One can see that the compositions are very different in two cases. 
In particular, a relevant hyperon (S~) fraction is only present in the MC. Non-uniform 
structure and suppression of hyperons in the mixed phase should have some implications 
on the neutrino transport and cooling of pulsars. 
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